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DEPENDENCE OF THE MAGNETIZATION 

OF AN ENSEMBLE OF SINGLE-DOMAIN PARTICLES 

ON THE MEASUREMENT TIME WITHIN VARIOUS EXPERIMENTAL 
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The effect of a measurement time duration on the parameters of magnetization curves for an 
ensemble of identical noninteracting single-domain particles with equally oriented axes under the 
uniaxial anisotropy has been specified for different experiment modes, in particular for the cases 
of relaxation measurements and the continuous sweep of a static magnetic field. The relation 
between a blocking temperature and experiment characteristics has been found for these modes. A 
recursion method to calculate the magnetization reversal curves for such an ensemble of particles is 
proposed. By comparing the results of calculations of the magnetic properties by the recursion and 
Monte-Carlo methods, an algorithm to establish the relation between the equivalent measurement 
time and such parameters of the Monte-Carlo method as the number of steps and the value of 
aperture is suggested. 
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I. INTRODUCTION 

The problems of the magnetism of nanoparticles have 
attracted the attention of scientists for many decades. 
More than a half-century ago, the English physicists E. 
Stoner and E. Wohlfarth developed a simple model of 
magnetization reversal for uniaxially anisotropic single- 
domain particles at T — K [1]. According to this model, 
during the reversal, all spins in each particle turn in such 
a way that they remain parallel to each other all the time, 
i.e. the absolute value of the magnetic moment for each 
particle remains constant, and only a mutual orientation 
of the magnetic moments of various particles changes. 
In such a case, the energy of the ensemble of particles 
depends only on one collective variable, for example, on 
a total magnetization vector. Within the frame of this 
model, all the particles constituting a sample are assumed 
to have the same shape, volume V, and orientation of the 
crystallographic anisotropy axes. The volume fraction / 
of these single-domain particles in a specimen is suffi- 
ciently small and, thus, the interparticlc interaction can 
be neglected. The objects, whose behavior corresponds 
to this model, are small (in order to satisfy the condi- 
tion for them to be single domains) magnetic particles 
(see, for example, [2]) placed in a nonmagnetic metallic 
or dielectric matrix. 

The density of magnetic energy U for a sample can be 
represented as the sum of the energy density of a mag- 
netic anisotropy, which includes the anisotropy caused 
by demagnetization fields (a shape-dependent anisotropy 
term) and that of the interaction of the magnetic mo- 
ment with an external field H. In the simplest case of 



a uniaxial anisotropy with regard for the first anisotropy 
constant only, the density of magnetic energy has the 
form 
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Here, Ui is the density of magnetic energy for the sepa- 
rate i-th particle, N is the number of magnetic particles 
in the sample volume, K is the first constant of the uni- 
axial anisotropy of a particle, m is its magnetization, 9i 
is the angle between the magnetic moment of the i-th 
particle and the magnetic field direction, and 6\ is the 
angle between the easy magnetization axis of a particle 
(for all the particles its direction is assumed to be identi- 
cal) and the magnetic field direction. For a sample which 
resides without an external magnetic field as long as pos- 
sible, the total magnetization M will turn to zero, since 
the numbers of particles with the magnetic moments ori- 
ented in parallel and antiparallel to the easy magnetiza- 
tion axis equal each other. On the contrary, the value of 
magnetization becomes finite upon the application of a 
magnetic field to the ensemble of particles. Hence, in the 
case where the observation duration of such the ensem- 
ble is infinitely long, its behavior will be characteristic of 
a paramagnet, in spite of the fact that the particles are 
ferromagnetic. 

Neel [3] and Brown [4] took the fact into considera- 
tion that a drastic removal of the magnetic field from the 
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ensemble of single-domain particles results in a time de- 
cay of the residual magnetization M(t) according to the 
exponential law: 



M(t) = M(0) exp - 



t 



(3) 



where M(0) is the magnetization value at the initial time 
moment and r is the relaxation time. The latter charac- 
terizes the thermally activated reversal of the direction 
of the magnetic moment of a separate particle between 
two possible minima of its potential energy. The hopping 
probability obeys the Arrhenius law which yields 
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where k is the Boltzmann constant, T is the tempera- 
ture, and E m i, E m 2, and E m b are the energies for two 
minima and a barrier between them, respectively. The 
expressions for these energy values contain the product 
of the particle volume V and the corresponding energy 
density Ui which is dependent on the relative orienta- 
tion of the particle magnetic moment and the field. As a 
result, the relaxation time r strongly depends on the par- 
ticle volume, temperature, and the value of the applied 
magnetic field. If the magnetic field is zero, the energy 
minima are equally deep. The application of a finite field 
makes the energy of one of them increase, while that of 
another one decrease. In a strong magnetic field, when 
the particle's Zeeman energy exceeds that of the uniaxial 
anisotropy, the higher-energy minimum disappears. For 
the magnetic particles under discussion, the typical val- 
ues of the prcexponcntial factor i>q in (|U) are between 
10 s and 10 10 s _1 . In the Arrhenius law, this factor 
is called "the attempt frequency". For estimations, its 
value can be assumed to equal the precession frequency 
for the magnetic moment in an effective magnetic field. 
Neel named these materials, which are the independent 
single-domain magnetic particles, superparamagnets and 
called their quasiparamagnetic behavior as superparam- 
agnetism. 

For the superparamagnets, the shape of magnetization 
curves strongly varies, depending on the duration t of the 
measurement process (measuring time). For each chosen 
value t, a blocking temperature Tf, can be introduced, 
which divides the whole region of temperatures into two 
ones with different magnetization behaviors. For one of 
them, the hopping, which occurs during the measuring 
time, of the particle magnetic moments between two en- 
ergy minima should necessarily be taken into account. 
But, for the second one, these effects are not essential 
and, thus, can be neglected. It is suitable to choose the 
temperature, at which the temperature-dependent relax- 
ation time r becomes equal to the measuring time t, as 
a blocking temperature. For T ^> T&, the measuring 
time i > t, and the magnetic moment of a particle has 
enough time to make multiple jumps between the energy 
minima. As a result, the populations of these minima 



do not differ from the equilibrium ones and the behav- 
ior of such a system of particles will be close to that of 
the ensemble of paramagnetic atoms, which is character- 
ized by the absence of magnetization hysteresis. In this 
case, the magnetization of the ensemble of particles can 
be described by formula 



M(H) — f m(H), 



(5) 



where the time-averaged particle magnetization, which is 
identical for all the particles, equals 
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The short-time measurements, for which i < r, corre- 
spond to T < Tft. In this case, there is no enough time for 
the transitions between the energy minima to occur, and 
the equilibrium populations for the states with different 
orientations of the magnetic moments of particles are not 
achieved during the measuring time. The system is in a 
metastable state and the curves of magnetization rever- 
sal display the hysteresis. The coercivity H c depends on 
the measuring time, anisotropy energy, and temperature. 
For the ensemble of uniaxial single-domain particles, Neel 
and Brown suggested a simple formula which connects 
these three parameters: 



H c (T) = H c (0)< 
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where t = v^ 1 , E = KV is the height of the energy bar- 
rier between the two minima at H = 0, and the exponent 
a is the parameter which depends on the orientation of a 
magnetic field relative to the easy axis of magnetization. 
For the case of the ensemble of uniaxial particles, whose 
easy axes are aligned along the magnetic field, a = 2.. 
For an arbitrary, but identical for all the particles, ori- 
entation of easy axes (with respect to the magnetic field 
direction), a is given as (see [5]) 

a = 0.86 + 1.14 ((cos 2/3 (0i) + sin 2 / 3 ^))-' 3 / 2 ) . (8) 

If the directions of the particles' easy axes are uniformly 
distributed over the space, a = 4/3. 

In the literature sources, for the case of a = 2, one can 
often find a representation of formula in the form 



H C {T) = H c (0) 1 - Jt/T™) 
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is the blocking temperature in the Neel-Brown approxi- 
mation. It is worth noting that, depending on the form 
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of the anisotropy energy definition in the above expres- 
sions, the anisotropy constant K can appear with a factor 
of 2. In spite of a relative simplicity, formulas - ([9]) 
successfully discribe experimental results. 

As was shown in work [6], formulas (O - © can 
be supplemented with the expression accounting for 
the dependence of the particle saturation magneti- 
zation on temperature. This will result in a de- 
viation of these expressions from the power depen- 
dence for the temperatures which are too close to the 
Curie temperature of ferromagnetic particles. The fur- 
ther improvements are reduced to the account of the 
distribution of particles over sizes or anisotropy val- 
ues or the account, in various approximations, of a 
magnetic dipole-dipole interaction between the parti- 
cles. 

However, the analytical calculation of the magnetiza- 
tion hysteresis curves for the systems under considera- 
tion meets serious difficulties even for a minimal number 
of independent parameters. At the same time, the power 
of modern computer systems makes it possible to carry 
out such calculations by numerical methods. One of the 
difficulties, which arise when one carries on the numer- 
ical calculations and tries to compare their results with 
experiment, is related to the correctness of the identifica- 
tion of a measuring time t, which appears in calculations, 
with a real duration of experiment. This implies that a 
relevant "protocol" of measurements should be taken into 
account. 

In modeling the properties of the ensembles of mag- 
netic nanoparticles, the method of Monte-Carlo (MC) 
[7-12] has gained a significant popularity. However, this 
method does not include the " real" measuring time. In- 
stead of it, the MC method contains such parameters as 
the number of mathematical iterations (MC steps) and 
the magnitude of angular aperture used to update the 
magnetization direction. At the same time, the liter- 
ature sources known to us do not contain an explicit 
relation between these parameters of the MC method 
and the equivalent measuring time, which corresponds 
to these calculations. In a few works (see, for exam- 
ple, [10, 13]), to correlate the number of MC steps with 
the measuring time, the results of MC calculations are 
compared with the data of actual magnetostatic mea- 
surements. The conclusions of these works are reduced 
to that the number of MC steps, being optimal from 
the viewpoint of the likelihood of the results obtained 
and the reasonable duration of calculations, corresponds 
to unrealistically short measuring times in real exper- 
iments, even for the calculations carried out with the 
use of high-performance modern computers. Though a 
number of papers devoted to this problem has been pub- 
lished for the last decade, the methods how to estab- 
lish the correspondence between the equivalent measur- 
ing time and MC simulation parameters remain ambigu- 
ous. 

It should be noted that the calculation, which would 
be able to account for the mode of carrying out the ex- 



periment, of the magnetization curves for the ensemble 
of single-domain particles has remained a problem, for 
a solution of which various approaches continue to be 
proposed (see, for example, [14, 15]). 

A method developed in this work for the modeling of 
the magnetization curves doesn't suffer from the above 
disadvantage. In what follows, we call it a recursion 
method (RM). In the literature sources, we haven't met 
any examples of the use of such a method. Its adequacy 
is grounded on the favorable outcome of the comparison 
of its results with those of both the MC simulations and 
basic formulas of the Neel-Brown model [3,4] described 
above. Basing on such comparison, we will be able to 
establish a specific relation between the parameters of 
MC simulation (the number of MC steps and the magni- 
tude of angular aperture) and the measuring time which 
corresponds to these parameters. The method we offer 
is suitable for the analysis of the magnetization curves 
for superparamagnetic systems consisting of uniaxially 
anisotropic particles and comprises the cases where the 
anisotropy axes are either parallel to each other or ran- 
domly oriented. There are no restrictions on its utiliza- 
tion to the modeling of the behavior of uniaxial systems 
with a nonzero second anisotropy constant and even the 
systems with a cubic anisotropy. However, to simplify the 
analysis, the modeling is carried out, in what follows, for 
the ensemble of uniaxial particles, whose axes are aligned 
in parallel to the magnetic field direction (0\ = 0) with 
regard for only the first anisotropy constant. 

II. DEPENDENCE OF BLOCKING TEMPERATURE 
ON MEASURING TIME FOR VARIOUS 
MEASUREMENT PROTOCOLS 

Consider an ensemble of identical noninteracting 
spherical single-domain particles, each of which has vol- 
ume V and is characterized by the uniaxial crystallo- 
graphic anisotropy. We assume that the easy axes of 
particles are aligned in parallel to the external magnetic 
field, i.e. we choose d\ = 0. We take only the first 
anisotropy constant into account. Let us rewrite formula 
@ in terms of dimensionless units by carrying out a divi- 
sion of its left and right parts by the anisotropy constant 
K 

U ied = -cos 2 (6») - 2hcos(0), (11) 

where h — Hm/ (2K) is the dimensionless magnetic field. 
Here and below, the index i in the notations of the energy 
density for a separate particle and the angle characteriz- 
ing the direction of its magnetic moment will be omitted. 
The parameter T rec i = kT/(KV) is used as a dimension- 
less temperature. Let us also give the definition of the 
dimensionless measuring time t rc d = tv$, where vq is the 
parameter which has a frequency dimension [see expres- 
sion (TJI for the probability of the thermally activated 
reversal of particle magnetic moments] and t is a real 
measuring time in seconds. 
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For the fields h £ (— 1, 1), the solutions of the equation 
dl gg d = give us the coordinates of two energy min- 
ima: 9 m i — n and m i = 0. The barrier between these 
minima is observed at 9^ — arccos(— h). The reduced en- 
ergies corresponding to these angles are E\ = (— 1 + 2/i), 
E 2 = (—l — 2h), and Et = h 2 . Substituting these quanti- 
ties into formula ((4J, the expression for the dimensionlcss 
relaxation time r rc d, which characterizes the thermally 
activated jumps between these minima, can be written 
as 



7"rcd 



As was noted above, the blocking temperature Tf, de- 
pends on the measuring time < rc( j- It is seen that the use 



of the condition t 
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17) for the determi- 



nation of the relaxation time leads to the relation 



T'b = l/ln(2t red ). 



(13) 



Here, TT is the blocking temperature in a zero magnetic 
field taken in the dimensionlcss form defined above. This 
expression doesn't coincide with that for the blocking 
temperature in the Neel-Brown approximation, T[, (a 
reduction of T b NB determined from expression (|10p to a 
dimensionless unit results in T b r » = l/m(£ re d))- Thus, the 
Neel-Brown approximation corresponds to neither the 
condition t = t(T = Tf,) (with expression for r in the 
form Q) nor this relation. 

Actually, formula © along with expression PH|) for 
T b NB can be obtained from the following considerations. 
In the case of high fields and low temperatures, i.e. when 
h/T led ^> 1, we can neglect the second exponential term 
in (|3]) (or in a dimensionless expression (fT2)) ). Let us 
take into account only the time, which is necessary for 
the thermally activated reversal of a particle magnetic 
moment from a metastable to the basic state, and ignore 
the backward jumps. In this case, 
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(14) 



A formal extrapolation, which is not strictly accurate, 
of this expression to the zero magnetic field and its equat- 
ing with a measuring time leads to a definition of the 
effective blocking temperature T b r „ in this approximation 

as 



T[» = 1/lny, 



(15) 



i.e. to the formula which was obtained by means of re- 
duction of (fT0|) to the dimensionless units. It is seen that 
77 (see (fT3"|) ) and TL (see ([T5]) ) are connected by a simple 
relation (T^)- 1 = (77) _1 - ln(2). 

The most important point in the approximation [3,4] is 
probably the expression for the temperature dependence 
of coercivity [see ([7]) and ©]• At low temperatures, the 
criterion h/T led ^> 1 can already be fulfilled for a coer- 
cive field [h — h c ) and, thus, approximation (|14p . along 



with the definition, which follows from it, of the blocking 
temperature, i.e. 77, becomes justified. 

Consider the question as to which type of experimen- 
tal data and measurement mode corresponds the Neel- 
Brown approximation in more details. Expression ([9]) 
with definition (JTUJ) [or (|12[) in dimensionless units] for a 
blocking temperature can be obtained proceeding from 
the assumption that the field which corresponds to coer- 
civity is the one, for which T re d(h = h c ) = t Ie ^. Then the 
relation 
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will be valid for the low temperature region of the coer- 
civity vs temperature dependence. After taking the log- 
arithm on both left and right parts of lfT6)) , we obtain the 
expression for the temperature dependence of coercivity: 



h e (T, 
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(17) 



Defining a blocking temperature as that, at which the 
extrapolation of a low temperature part of the tempera- 
ture dependence of coercivity with either formula |9]) or 
(|17|) reaches zero value for a preset measuring time, we 
can write the expression for the temperature dependence 
of coercivity as 



^c(Trcd) — 1 — Y^Tred/lb*, 



(18) 



where 77* is determined from formula (|15p . 

However, one should keep in mind that, to derive Eq. 
(fT8)l . we had to use assumption |[T6)) which is not com- 
pletely accurate in the strict sense. Moreover, it should 
be remembered that, to obtain (|15p . we used one more 
approximation which consisted in the extrapolation of 
the low-temperature high-field part of the r* cd vs h de- 
pendence to the zero field. In fact, the measuring time, 
which goes into these equations, can be interpreted as 
the time, during which the system relaxes from the mag- 
netosaturated state after the instantaneous switching-on 
of the given field. Such a definition of the measuring 
time is related to the relaxation measurements. For the 
relaxation experiments under consideration, the time de- 
pendence of the magnetization is given as 

M(h, t re d) = M equ (/i) + (M sign(/i satur )- 



M equ (h)) exp 



-(/i,T red ) 



(19) 



Here, M eqn (h) is the equilibrium magnetization in a field 
h, which for the system under study is determined as 
M cqu (h) = M tanh(2/i/r rcd ), sign(/i satur ) is the sign of 
the saturation field, and Mq is the saturation magnetiza- 
tion. At the same time, the coercivity is determined from 
the condition M(h c , t Te d) = which differs from the con- 
dition t Te d = T rcd (h Cl T red ). The latter, perhaps, might be 
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used as an approximate condition. At low temperatures 
and sufficiently high coercivities (when h c /T Ie d 3> 1, 
which is a criterion of the applicability of approximation 
tH}), it is either a relation exp [— t Te d/T Te d(h c , T rc( j)] ~ 
1/2 or £ rc d = T ro d(/i C 7 Trcd) hi(2) that would better satisfy 
the requirement M(h c , irod) = 0. Thus, we will use them 
instead of (|16p . The utilization of such an approximation 
makes it possible to obtain the equation for the temper- 
ature dependence of coercivity. This equation has the 
same form as formula ©, in which the effective block- 
ing temperature T£, is substituted by Tgi*. The latter 
coincides with neither ([TUf nor fTS}; it equals 



b " ln[i rcd /ln(2)] ln(t rcd ) - ln[ln(2)] ' V ; 

It is seen that (7^ J" 1 - (I^)" 1 - ln[ln(2)] =7J - 
ln(2) -ln[ln(2)]. 

To correctly determine T& and obtain the expression 
for the coercivity vs T dependence, one should not only 
find the measuring time, but also take the measurement 
protocol for the magnetic characteristics of a system of 
particles into account. 

Consider the relaxation experiments where the satu- 
rating magnetic field, which acts on the system, is in- 
stantaneously (or rapidly enough in real conditions so 
that t < t) transformed into some other field h, after 
which the system relaxes with time. In this case, let the 
measurements of the magnetization be continuously car- 
ried out during a time interval t Ie d- Consider the ways 
of the determination of both T[ [see (|13p] and the ef- 
fective blocking temperature T 6 r » [see (|18p] to make it 
possible to specify the temperature dependence of coer- 
civity. It is pertinent to take a temperature at which 
T le d(h = 0) = t rc d [with regard for (fT2"|) for r ro d] as T b r . 
If this condition is fulfilled, the zero-field magnetization 
(the so-called remanent magnetization) has to be e times 
less than the saturation magnetization (here and below, 
the number e means the base of the natural logarithm). 
To measure the coercivity according to this method and 
to determine the blocking temperature from the h c (T TC d) 
dependences, it is necessary at each measurement tem- 
perature and at a certain field which has the opposite sign 
relative to the initial saturation field, to record a time 
point, at which the magnetic relaxation curve crosses the 
zero value. The corresponding field will equal the coer- 
civity at the given values of temperature and measuring 
time. At low temperatures, it will correspond to the so- 
lution of the equation 

, r f (1 - h c (T Icd ,t Tcd )) 2 \ 
t lcd = [cxp I — I + 

+ ex P Mi + MWrcd)) 2 ^-!^ (21) 
\ T T ed J 

Then, for each value of t Te d, one should plot the depen- 
dence h c (T re d)\t TCd and find (t ie d) from the approx- 
imation of its low temperature part by formula (|18p . 
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FIG. 1 Temperature dependences of coercivity calculated for 
different measuring times: f ro( j = 4, 10, 10 3 , and 10 7 for the 
pairs of curves 1, 2, 3, and 4, respectively. Dashed and solid 
lines represent the solutions of Eq. (|21|l and M(h c ,t le d) = 0, 
respectively. 



in which XJ^ substitutes for TjJ"„. A solution of Eq. 
([2"Tjl . which describes the temperature dependence of co- 
ercivity, can be found only numerically. As was noted 
above, it gives a correct criterion for the determination of 
h c (T re d, ired) only for the temperatures lower than T^. 
To find the temperature dependence of h c over the whole 
range of temperatures, it is necessary to solve the equa- 
tion M(h c , i ro d) = with the use of expressions (TT2"|) 
and (fT9|) for r rc d and M(h,t re d), respectively. It should 
be noted that the h c (T Te d) dependence found in such a 
way doesn't exhibit a sharp break by turning into zero 
at the blocking temperature. On the contrary, it dimin- 
ishes smoothly over a certain temperature range above 
the blocking temperature. 

Figure 1 shows the temperature dependences of co- 
ercivity for various measuring times t Te d calculated by 
means of the numerical solution of Eq. (f2"Tj) and 
the equation M(h c , t le d) = with the components 
described above. The quantities y/T TC d/T^ (i re d) = 
\/2red • {hi(trcd) — In [l n (2)]} are taken as the units of the 
abscissa axis. 

As is seen from the figure, all the curves coincide 
with each other at low temperatures. Extrapolation of 
the low temperature branches to higher x values gives 
y/ Tl-od / ~ ( t TO d ) = 1- However, the real curves do not 
follow the extrapolated one. On the contrary, they di- 
verge to different sides in the vicinity of the effective 
blocking temperature. The curves obtained from Eq. 
([2"Tj) turn into zero below , while those obtained from 
the condition M(h c ,t re d) — diminish smoothly above 
Tfc**- For very great values of i re d (£red > 10 5 ), the curves 
remain linear practically to T^. It should be noted that, 
taking i/g w 10 s ~ 10 10 s _1 into account, only such val- 
ues of t ro d are characteristic of real measurements. The 
shorter the measuring time, the lower is the tempera- 
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ture, at which the curves start to deviate from the linear 
law which corresponds to formula (|18p where the effec- 
tive blocking temperature equals T^. At the same time, 
both the temperature smearing of a sharp transition in 
the vicinity of h c = and the transition retention to the 
temperatures above the blocking temperature disguises a 
deviation of the blocking temperature from Tjk* at short 
(almost unlikely in practice) measuring times. For this 
reason, if one employs the relaxation method to measure 
a coercivity, the deviation from formula (|18p . in which 



T£.„ serves as the effective blocking temperature, will 
appear only in the form of the aforementioned smear- 
ing. The temperature smearing of the transition, which 
is observed as h c tends to zero, seems to be natural, since 
the hysteresis is a manifestation of the metastability at 
a finite measuring time. Such a hysteresis will occur, to 
a greater or lesser extent, at any finite temperatures. It 
is worth noting that the coercivity, even being negligibly 
small at T > doesn't turn into zero. 

In principle, the relaxation measurement protocol con- 
sidered above is used in practice. However, the proto- 
col of the continuous sweep of a magnetic field (CSMF) 
with certain rate is more often used for magnetostatic 
measurements. In the course of its implementation, the 
relaxation of the magnetization to its equilibrium value 
occurs in a magnetic field, which continuously changes. 
The case where the sweeping rate is infinitely small cor- 
responds to the infinitely great measuring time. In this 
case, the system is in an equilibrium state, and such a 
case corresponds to T\, — > 0. For the regions where super- 
paramagntism becomes clearly apparent and Xf, ^ 0, the 
sweeping rate becomes comparable with the relaxation 
time t. Under these conditions, the concept of a mea- 
suring time should be made more specific and related to 
the experiment conditions and the blocking temperature 
definition. 

In the case of the CSMF studies, the simplest and nat- 
ural way to analyze and describe the system properties 
may be the analysis of a hypothetical protocol of mea- 
surements, in which the whole range of the field sweep is 
divided into equal intervals. In this case, the magnetic 
field sweep can be regarded as a series of jumps, each of 
which being characterized by a specific waiting time t w 
after the previous jump. In this case, the magnetization 
for each field point will relax during t w from the value 
at the previous point. Such a method can be called as 
recursive, since, in order to describe the magnetization 
relaxation at the n-th field point, one should reconstruct 
a successive series of magnetization relaxations for all 
n — 1 previous points. In the limit where the interval 
between successive points tends to zero, we obtain the 
CSMF protocol. It is appropriate to define the measur- 
ing time as a sweep time for the unit field interval (taken 
in dimensionless units), i.e. to define i rc d as the quantity 
which is reciprocal to the sweep rate averaged over the 
whole intervals. 

Figure 2 shows the dependences of the blocking tem- 
peratures defined in different ways on 1/ ln(t rc( j): T£ was 




0,18 



FIG. 2 The Tb_ scan (£ rc d) dependences calculated by the RM 
technique for various numbers of points per unit field interval 
(circles). The top solid line corresponds to the limiting case 
n — > oo. The bottom line shows the T 6 r (£ rc d) dependence 
calculated by formula (|13[) 



obtained from expression (JT3J) , while T£_ scan was calcu- 
lated numerically by means of the RM for different values 
of n, where n is the number of intervals, into which the 
unit field interval Ah = 1 was divided. In calculations, 
the ?b_ scan was regarded as a value of the reduced tem- 
perature, at which the remanent magnetization was e 
times less than the saturation magnetization. 

It is seen that all dependences, which are calculated 
with the use of the RM, lie between two curves, one of 
which is the curve obtained from (fT5)) . while the other is 
some limiting curve, to which the calculated results tend 
when n goes to infinity. This curve corresponds to the 
CSMF case. It is clear that the solution corresponding 
to n — 1, i.e. when one point falls on a unit interval, 
coincides with that obtained from formula (|13[) . One can 
see that the results start to diverge from the limiting 
solution already for n ~ 200. 

The dependences h c (T rc a) were also calculated for the 
case n — > oo. Extrapolation of their low temperature re- 
gions to the temperatures where h c — > resulted in the 
same dependences of T£_ scan (£ ro d) as those which were 
determined from the expression for the remanent magne- 
tization. For all t re d values, except for the shortest ones 
(i rc d < 5), the curves /i c (T re d) correspond to Eq. |T5J) , in 
which T£_ scaa (t re d) substitutes for T b r v 

The analysis of the data in Fig. 2 shows that the CSMF 
case is the limiting case of relaxation measurements when 
the measuring time is reduced. In fact, in the former case, 
the system resides in the saturated state for a consider- 
able part of the sweep time which falls on a unit field 
interval. The effective sweep time from the saturating 
field to the zero (or coercive) 

one turns out to be far shortest than t le d defined above for 
the CSMF protocol. For this protocol, one can introduce 
the effective measuring time t^ d , at which T£_ scan (£ re d) 
in the CSMF case will coincide with T£(t°^ d ) in the case 
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of relaxation measurements, i.e. it will equal l/ln(2^ d ). 
It is appropriate to regard such time as the sweep pe- 
riod for the field, at which the equilibrium magnetization 
changes from zero to a certain value clMq. Then, from 
the fitting of the limiting curve in Fig. 2, we should find 
the optimal value of a (its expected value is about 2). 
Such an approach leads to the equation t°^ d = i rod arctg 
h(a)/ [21n(2^f d )] , whose solut ion is 



ired(^red) a ) 



i re d arctanh (a) 
2W [t rc( j arctanh (a)] 



(22) 



Here, W(x) is the so-called Lambert PF-function which 
is a reciprocal function to x = W(x) exp(W(x)). This 
function was introduced into mathematical physics rela- 
tively recently [18] . In a number of cases, for example in a 
popular program "Mathematica" developed by "Wolfram 
Research" company, it is denoted as "ProductLog" . The 
examples of the solution of various tasks of mathematical 
physics with the use of this function are presented in [19]. 

The calculations we carried out have shown that the 
expression 



TT 



1/ In [2C(*rcd, a = 0.45167)] (23) 



approximates the limiting curve of Fig. 2 with a high 
precision. Noticeable deviations are observed only for 
the smallest i rcd values (t Te & < 5), which are of minor 
importance from the practical point of view. Thus, this 
approximation requires only one fitting parameter a. 

Taking into account that the W(x) function is not 
widely used in the scientific literature, we also found the 
approximation of the limiting curve for T^_ scan (t ro d) by 

a power series of [ln(2t rcd )] : 



1 + a + 



rpr 

b— scan 



Ol 



{tied) 



1 



In (2t red ) 



«3 



ln(2t rcd ) [ln(2i rod )] 2 [ln(2t red )]3 

(24) 

with four fitting parameters: — 0.01, a\ — 5.197, a 2 = 
—3.581, and 03 = —1.602. Finally, we note that, in real 
CSMF experiments, £ rcd is usually between 10 10 and 10 14 . 
For this reason, the discrepancies in the determination 
of a blocking temperature are almost unnoticeable for 
different measurement methods. 



III. RECURSION METHOD 

Before turning to a description of the calculation 
method for magnetic reversal curves for different mea- 
suring times, it should be noted that its implementation, 
as also in the case of MC modeling, requires some efforts 
in the programming of a calculation procedure. We will 
use the same dimensionless parameters as in Section 2: 
h = Hm/(2K) is the magnetic field, t rod = tv$ is the 



measuring time, T red = kT/(KV) is the temperature, 
n is the number of points per unit field interval, and 
m red = fm/(fm s ) — M red - is average magnetization for 
measuring time, which normalized on a saturation mag- 
netization m s . 

The method is based on the following prerequisites. 

1. The hysteresis, which is a result of the metastabil- 
ity, will become apparent at a finite measuring time only 
if the state of a system is characterized by two minima in 
the dependence of its energy on the orientation of a parti- 
cle magnetic moment. The hysteresis originates from the 
metastability with regard to thermally activated jumps 
over a potential barrier. In its turn, the metastability 
appears as a result of the finiteness of a measuring time. 

2. In real magnetostatic measurements, the ratio of a 
barrier energy at H = to a thermal energy, at which 
the deblocking of a magnetic moment occurs, is near 25. 
That is why we assume that, even at temperatures higher 
than the blocking temperature, the orientation of a mag- 
netic moment will be localized in one of the minima, 
rather than smeared by temperature over a wide range 
of angles 9. 

3. To determine the magnetization of a system, the 
concept of potential well (minimum) populations is intro- 
duced. This concept is based on the distribution statis- 
tics of magnetic moment directions in an infinitely great 
ensemble of identical and equally oriented particles. 

4. In dimensionless units, the magnetization of the 
system depends only on the coordinates of minima 
(6 m i,6 m 2) and their populations (JVijiVjj): M red = 
N lC os(9 ml ) + N 2 cos{9 m2 ). 

The limits of applicability of this method will be dis- 
cussed later on; we will concentrate now only on the cal- 
culation procedure. 

We reduce formula @ for the density of energy of a 
separate particle to that in the dimensionless units (the 
index i, which refers to the number of a particle, will be 
omitted again, as was done in (jXTJ) ) : 



U r = - cos 2 (0 - 0i ) - 2/i cos(0) . 



(25) 



Consider the energy profile U r in the phase space 9 S 
(— 7r, 7r] as a function of the magnetic field h. The energy 
minimum will correspond to the equilibrium orientation 
of the magnetic moment. The number of extrema can be 
found by solving the equation dll r /d9 = 0. Substituting 
the roots of this equation into the expression d 2 U r /d9 2 , 
we can find out if a given root corresponds to the energy 
minimum or maximum. 

The method developed consists of a few successive 
steps. At first, the system is assumed to reside in a sat- 
urating magnetic field. Then we sequentially change the 
field to smaller values and calculate the magnetization, 
to which the system will come during the time interval 
which is equal to the period of the system residence at 
a certain field point h^. The relaxation time is deter- 
mined by the form of the potential in a field . Accord- 
ing to this method, the calculations should start from a 
negative field, which is high enough so that the energy 
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displays only one minimum, and finish at a sufficiently 
high positive field, for which the energy again displays 
only one minimum. Thus, the calculation procedure can 
tentatively be divided into three stages. 

The first stage is applicable when the field is negative 
and the state of the system is characterized by only one 
energy minimum. This stage includes: 

1. Determination of the total number of extrema in 
the range 9 G (— 7T, tt] and their separation into minima 
and maxima. 

2. If the state of the system is characterized by only 
one energy minimum, then iVi — 1 and N 2 — 0, and the 
magnetization equals M ro d = Ni cos(# m i) + N 2 cos(9 m2 ). 
If there are two minima, we should go on to the second 
stage. 

3. Recording the magnetization for a given field point 
hk, changing the field to hk+i — hk + A, and going on to 
item 1. 

The second stage is applicable when the state of the 
system is characterized by two energy minima and in- 
cludes: 

4. Determination of the total number of extrema in 
the range 9 G (— tt, tt] and their separation into minima 
and maxima. 

5. If the state of the system is characterized by a single 
energy minimum, then we should go on to the third stage. 
If there are two minima, a temporal evolution of N\ and 
N 2 should be considered. This includes: 

a) determination of the equilibrium populations Niaa 
and A/200 > i-e. the values of Ni and N 2 when the time 
interval is infinite: 

^100 = — j— 77; „ \ ,„ — t> N 2oo = l-Ni 00 , 

1 + cxp [-(E 2 -E 1 )/T ied ] 

(26) 

where E\ — U T (9 ml ) and E 2 = U r (9 m2 ) are the energy 
values for the first and second minima, respectively; 

b) tracing the relaxation of N\ and N 2 for the time 
interval ii oc = t rec j/n: 

N u = N loo + (Ni - N loo ) exp {-t ioc /r led ) , (27) 

1 

Tr ° d ~ cxp [-(E b - Ei)/T red ] + exp [-(E b - E 2 )/T red ] ' 

(28) 

c) recording the new values of Ni and N 2 : 

N X =N U , N 2 = 1-Nu. 

6. Calculation of the magnetization Af rc( j = 
Ni coa(9 ml ) + N 2 cos(9 m2 ). 

7. Fixation of the magnetization for a given field point 
hi, changing the field to hi+i = hi + A, and going on to 
item 4. 

The third stage is applicable when the field is positive, 
and there is no second minimum. This stage includes: 

8. Determination of the total number of extrema in the 
range 9 S (—tt, tt] and separation of them into minima 
and maxima. 



9. If the state of the system is characterized by only 
one energy minimum, then iVi = and N 2 = 1 and the 
magnetization equals M ro d = Ni cos(# m i) +-/V2 cos(0 m 2). 

10. Fixation of the magnetization for a given field point 
h m , changing the field to h m +i — h rn + A, and going on 
to item 8. 

The greater n, the more minutely the magnetization 
curve will be calculated. At the same time, the real time 
interval spent on the measurement of the complete hys- 
teresis loop ({ — h, h}, {h, —h}) will equal t — 2t rc dA/i/^o 
s, where Ah is the interval (taken in dimensionless units) 
of the complete sweep of the magnetic field. If one 
does not need to model the partial hysteresis loops, it 
is enough to carry out the calculations for the interval 
{ — h,h}, because the second interval {h,—h} will be 
symmetric relative to the coordinate origin (0, 0). To ob- 
tain the magnetization curves for an ensemble which is 
characterized by a distribution of some particles' param- 
eters, it is necessary to divide the distribution function 
into sufficiently small intervals and, having calculated the 
separate curves for each of the intervals (i.e. for the 
average values of a parameter in this interval), to sum 
them. 

The inset in Fig. 3 shows the results of the mod- 
eling of magnetization curves for the ensemble of uni- 
axial particles, whose easy axes are aligned in parallel 
to the magnetic field. The calculations were made for 
^rcd = 1-25 x 10 6 and n = 125. The temperature T re d was 
varied from to 0.02. The temperature dependence of 
coercivity (squares) is well described by formula (0 and 
has a square-root character (a = 2) for the given orienta- 
tion of the particle easy axes. The plot of this curve in the 
form /i c (yT re d) (triangles) confirms the latter fact. Such 
a dependence is a straight line and its extrapolation to 
the intersection with the ordinate axis unambiguously de- 
termines the blocking temperature T£ ed = (0.3 ± 0.005) 2 
= 0.09 ± 0.003. At the same time, the calculation of 
T b'- S can(^cd = 1-25 x 10 6 ) according to formulas {23]) 
and (|24j) in the case n — > 00 (or, at least, n > 200) gives 
0.0915. Thus, the value of T b rod = 0.09 obtained from 
Fig. 3 well agrees with that found from (|23|) and |24|) . 

The temperature dependence of the remanent magne- 
tization M r normalized to the saturation magnetization 
exhibits a sharp rise (see Fig. 3, circles) at temperatures 
where h c becomes noticeable. 

As was expected, the value T ro d = T h rcd = 0.09 cor- 
responds to a decrease in the remanence by a factor of 
2.718. This means that the plots of the temperature de- 
pendences of the coercivity and the remanent magneti- 
zation normalized to the blocking temperature (see ex- 
pression (|23|) or (|24p) should coincide. Figure 4, a shows 
the results of the RM modeling obtained on the same 
ensemble of particles, but for different values of the mea- 
suring time. These data replotted in the coordinates 
where the abscissa axis is scaled by T£_ scan (£ re d) are 
presented in Fig. 4,&. It is seen that almost all the 
H c [T , re d/Tb_ scan (t re d)] curves overlap one another. Re- 
plotting these curves in the coordinates with a "rooted 
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FIG. 3 Results of the modeling of the h c (T re a) (solid squares) 
and M r (T re d) (open circles) dependences for t rc d = 1.25 x 10 6 . 
The h c (\/T TC d) dependence is shown by open triangles. The 
calculations are carried out for n — 125. The inset shows the 
magnetization curves which correspond to these dependences, 
for T rcd = 0.2, 0.1 (anhysteretic curves) and 0.07, 0.03, 0.01, 
and 0.00 (with a hysteresis which consecutively widens) 



abscissa" shows that, at temperatures far lower than the 
blocking temperature, they are well described by formula 
©■ 

A characteristic feature of the M r [T re d/7b_ scan (tred)] 
dependence is a broadening of the front of the M r growth 
with decrease in i xe d. All the curves intersect one an- 
other at a point, whose ordinate equals 0.368, i.e. 1/e. 
Thus, the simulated magnetization curves, as well as the 
H c (T le d) and M r (T re d) dependences obtained from them, 
are in compliance with the regularities described in Sec- 
tion 2, which determine their behavior over a wide range 
of measuring time values. 

The RM calculations of the magnetization curves were 
also carried out for the ensembles of particles with 3D- 
and 2D-distributions of easy axes. The results obtained 
well agree with the data calculated by the MC method 
[1,6,7]. The fact that the time of computer calculations 
is much less within the RM than that within the MC 
method gives us an additional argument in favor of the 
method we have proposed here. 

It is worth noting that the method developed contains 
a series of approximations, which can lead to some in- 
accuracy of the results obtained. The most important 
approximation is related to a failure to consider the ther- 
mal fluctuations of a magnetic moment in the vicinity of 
an energy minimum . This, in turn, gives rise to an in- 
accuracy in the determination of the magnetization at 
sufficiently high temperatures (or for very short values of 
the measuring time (£ re d < 5)). It should be stressed once 
more that, in the first place, the RM is used for the cal- 
culations of the ithysteresis loops of the magnetization. 
To study the magnetization curves itabove itthe blocking 
temperature, it is enough to utilize formula ([5]). With the 



use of this formula, one can also estimate the measure of 
inaccuracy for the magnetization calculated by the RM 
at temperatures higher than the blocking temperature. It 
is obvious that the higher the temperature, the greater 
is the inaccuracy. To check the role of this factor, we 
calculated the magnetization curves for t Ye a = 1-25 x 10 4 
and T rcd > T fc rcd (namely, for T rod = 0.2) by the RM and 
formula ©. The maximal error in the determination of 
the magnetization did not exceed 1.5%. The measuring 
times, which are shorter than the above value, are hardly 
possible in practice. 



IV. MONTE-CARLO METHOD 

For the modeling by Monte-Carlo method, the stan- 
dard algorithm suggested by Metropolis et al. [16] was 
used. It is known that, for a sufficiently great number 
of steps Nmc, such an algorithm leads to the Boltzmann 
distribution. This means that a system comes to the 
thermodynamic equilibrium and thus, no metastability 
and, respectively, no hysteresis will be observed unless 
we introduce some special tricks. In a general case, for a 
great number of MC steps, the results will tend to those 
which can be obtained with the use of formula ([5]). To 
"catch" the metastability in the process of magnetiza- 
tion reversal, it is necessary to use a finite number of 
MC steps and restrict the generation of a trial random 
orientation of the magnetic moment in the vicinity of the 
current orientation by a certain not great aperture AO, 
instead of the generation over the whole phase space. 
This trick is one of the standard MC techniques to model 
the hysteresis loops [7, 10, 12]. The procedure of mod- 
eling is divided into two stages: a thermalization of the 
system and a magnetization reversal process itself. 

The system thermalization is regarded as the proce- 
dure consisting from tens to hundreds of thousands of 
MC steps at high negative fields and sufficiently high 
temperatures. This procedure is aimed at bringing the 
system to a state which is equivalent to the thermal equi- 
librium. To model the magnetization reversal process, 
one should give an increment in the magnetic field and 
perform a given number of MC steps for each field point. 
Such algorithm is widely used, and we won't describe it 
in details (see its description, for example, in [10,12,17]). 
Here, we only note that, in our calculations, we set the 
aperture value AO = 6° (the role of the aperture will 
be discussed below). The thermalization procedure was 
carried out for T rod = 0.5, h = -10, and A M c = 10000. 
To compare the efficiency of the MC method and RM, 
we performed a series of calculations of the magnetiza- 
tion curves for various T rc d by both methods. To de- 
termine the blocking temperature, we utilized the ex- 
trapolation of a low-temperature region of the /i c (^/T ro d) 
dependence, which is linear in these coordinates, to its in- 
tersection with the abscissa axis. In both cases, 300 field 
points fell on one magnetization curve (—2 < h < 2), 
which meant that n was equal to 75. Then we carried 
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FIG. 4 (a) Families of the H c (T Ic a) (solid squares) and 
M r (T Ic d) (open circles) curves calculated for the ensemble 
of particles, whose easy axes are aligned along the magnetic 
field, for t rcd = 2.5 x 10 12 , 2.5 x 10 10 , 2.5 x 10 8 , 2.5 x 10 6 , 
2.5 x 10 4 , and 2.5 x 10 3 (n = 250). The curves are shifted 
towards the right side as t rc a decreases. (6) The same curves 
replotted in the coordinates where the abscissa axis is scaled 

by ^b-scan (tied) 

out the calculations of the H c (T re d) and M r (T re d) depen- 
dences according to the MC method with Nmc = 5 x 10 6 , 
a measuring time £ re d for the RM procedure was fit- 
ted in such a way that the dependences obtained agreed 
as much as possible with the results of the MC model- 
ing. 

Figure 5 presents the H c (T le d) and M r (T rc< \) depen- 
dences obtained by the MC method (squares and cir- 
cles) and RM (solid lines) for the ensemble of particles, 
whose easy axes are aligned along the magnetic field. It 
is seen from the figure that, in spite of the disadvantages 
of the RM, which were formulated in Section 3, both 
the methods give almost identical dependences. Since 
the duration of the MC calculations was sufficiently long 
(it took 46 min to calculate one magnetization curve), 
only a small number of particles (5 particles in the en- 
semble) was taken for calculations and this resulted in 
a noticeable data scattering. On the contrary, it took 
no more than a minute to make the RM calculations for 
one magnetization curve. In the case under considera- 
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FIG. 5 Dependences H c (T Te d) (solid squares) and M r (T re d) 
(open circles) obtained by the MC modeling (JVmc = 5x 10 6 , 
iloe = 8 x 10 3 , n = 75) and the RM (solid lines) 



tion, the number of MC steps at each field point was 
5 x 10 6 , while the RM with a fitting procedure described 
above gave t\ oc — 8 x 10 3 , which means that, for n = 75, 
^rcd = nt\oc = 6 x 10 5 . Since the parameter vq is of the 
order of 10 8 — 10 11 , the actual measuring time for such 
a kind of the MC procedure corresponds to 10~ 4 — 10~ 7 
s. This is, of course, an extremely small time interval. 
At the same time, it is pertinent to specify the 
interrelation between the number of MC steps and the 
aperture, on the one hand, and the real time interval, 



V. RELATIONSHIP BETWEEN THE PARAMETERS OF 
MC MODELING AND REAL MEASURING TIME 

It is appropriate to assume that, for the MC procedure, 
the equivalent real measuring time i re d should be propor- 
tional to tiTVmo For this reason, the results obtained in 
Section 2 were used to search for such a dependence. At 
first, the calculations of the blocking temperature T h red 
were carried out by both methods for different measur- 
ing times (different numbers of MC steps). 

Figure 6 shows the dependences obtained along with a 
theoretical curve which corresponds to formulas (f2"3")) and 
To make a comparison of the results to be more 
covenient, all the dependences were represented in the 
form T h red [ln(i sys /i fit )], where t sys is the effective mea- 
suring time characteristic of each method (for the RM 
and the theoretical dependence, t sys = t Ie< i; for the MC 
procedure, t sys = uNmc), and £fit is the fitting parameter 
specific for each procedure. It turned out that, for the 
RM procedure and the theoretical dependence, igt = 1- 
It was also confirmed that, for the MC procedure, firstly, 
the equivalent measuring time is actually proportional to 
nNMC, and, secondly, the coefficient of proportionality is 
about 450. It is seen that the points calculated by both 
methods agree well with the theoretical curve. Based on 
this, the relation between the number of MC steps and 
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FIG. 6 Dependences of T£ ed on the effective measuring time 
i S ys (the number of MC steps) for the MC method (triangles) 
and the RM (circles). The solid curve corresponds to the 
calculations by formulas (|23p and (|24|l 
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FIG. 7 The T£ ed {nM MC ) dependences calculated by MC 
method for various values of aperture AO = 1.5, 3, 6, 12, 24, 
and 48° (circles). Solid lines are the fitting of the obtained 
results with the function Tb-scan^A^Mc/ifit) 



the equivalent measuring 
time was obtained in the form 



t 



MC 
red 



niV MC /450. 



(29) 



One can estimate the number of MC steps which cor- 
responds to the measuring time tv§ = 10 10 . It equals ap- 
proximately 5 x 10 12 . Taking into account that one MC 
step contains a few tens of operations and that a mod- 
eling is performed for an ensemble of particles, the nec- 
essary computational resources exceed 10 14 operations. 
Thus, it is concluded that the computational capability 
of modern computers is not sufficient to carry out the 
real time modeling. It is obvious that, for some calcula- 
tions, it is more important to deal with a great number of 
particles in an ensemble and less important whether the 
number of MC steps is great or small. However, to sim- 
ulate, for example, the ZFC/FC (zero field cooling/field 
cooling) procedure, the value of the measuring time is 
likely to be decisive. 

Let us discuss the origin of the parameter tat for the 
MC procedure. The only parameter which we have not 
yet varied is the generation aperture AO. It was noted 
that the MC calculations above were carried out under 
the condition that AO = 6°. Let us start from the as- 
sumption that it is the parameter AO that determines 
the ifit value. To make sure of this, we performed the 
MC calculations of T b red (A/Mc) dependences for various 
values of aperture (A9 = 1.5 48°). 

Figure 7 presents the results of this modeling. As is 
seen from the figure, for the calculations with a small 
number of MC steps, X£ ed strongly depends on Nmc- As 
Amc increases, the aperture value to a lesser extent af- 
fects T h red . At the same time, a family of T fc rcd [ln(nN MC )} 
dependences can be well approximated by a function 
7b-scan(ft-^McAfit)- f* follows from general considera- 
tions that tfa(AO) should be dependent on [/,., its deriva- 
tives, or its integrals. However, we have not succeed in 



finding the precise analytic dependence. We can only 
note that, for AO < 45°, the admissible expression re- 
sulted from the fitting of the data in Fig. 7 is 



tea « 4.9/ [sin 2 (A0)] 



(30) 



Thus, we have obtained the final expressions which an- 
swer a question about which measuring time corresponds 
to the MC calculations: 



.MC 

''red 



\A0) 



4.9 



nN MC , 



T, 



MC 



Nmc . 

n sin (AO) 

4.9 V ' 



(31) 



(32) 



VI. CONCLUSIONS 

In this work, the recursion method has been developed 
for the calculations of the magnetic properties of the en- 
semble of single-domain particles. Its applicability to the 
system of oriented particles with a uniaxial anisotropy 
is demonstrated. There are no hindrances to apply such 
a procedure to the case of a cubic anisotropy, introduce 
the distribution function for a certain parameter, make 
the anisotropy constant temperature-dependent, or even 
introduce the dipole-dipole interaction between the par- 
ticles of an ensemble. It is also not difficult to model the 
ZFC/FC procedure. In our opinion, the RM results will 
far better reflect the real experiments, than the results 
of the MC modeling. 

The relation, which correlates the magnetic parame- 
ters of the ensemble of uniaxially anisotropic magnetic 
particles with a measuring time for these properties in 
various experimental procedures, is obtained. 
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It is shown that, depending on a kind of experiment, 
the relationship between the blocking temperature and 
the measuring time has somewhat different forms. 

The calculations of the magnetization curves for the 
ensemble of uniaxial single-domain particles are carried 
out for different measuring times. The similar calcula- 
tions performed by the Monte-Carlo technique confirm 
the adequacy of the method developed here. The latter 
method requires far less computational resources 
in comparison with the modeling of an analogous task by 
the Monte-Carlo method. 

With the use of the new method, we succeeded in the 
establishment of the empirical dependence between the 
number of MC steps and the generation aperture of a ran- 
dom direction of the magnetic moment, on the one hand, 
and the measuring time, which corresponds to these pa- 
rameters, on the other hand. It is shown that the param- 
eters, which are usually used for the MC modeling of the 
behavior of ensembles of magnetic particles, correspond 
to unlikely short values of the measuring time. 

Translated from Ukrainian by A.I. Tovstolytkin 
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